Missense3D-TM: Predicting the Effect of Missense Variants in Helical Transmembrane Protein Regions Using 3D Protein Structures.

Variant effect predictors assess if a substitution is pathogenic or benign. Most predictors, including those that are structure-based, are designed for globular proteins in aqueous environments and do not consider that the variant residue is located within the membrane. We report Missense3D-TM that provides a structure-based assessment of the impact of a missense variant located within a membrane. On a dataset of 2,078 pathogenic and 1,060 benign variants, spanning 711 proteins from 706 structures, Missense3D-TM achieved an accuracy of 66%, Mathews correlation coefficient of 0.37, sensitivity of 58% and specificity of 81%. Missense3D-TM performed similarly to mCSM-membrane: accuracy 66% vs 61% (p = 0.02) on an unbalanced test set and 70% vs 67% (p = 0.20) on a balanced test set. The Missense3D-TM website provides an analysis of the structural effects of the variant along with its predicted position within the membrane. The web server is available at http://missense3d.bc.ic.ac.uk/.

GWAS of random glucose in 476,326 individuals provide insights into diabetes pathophysiology, complications and treatment stratification.

Conventional measurements of fasting and postprandial blood glucose levels investigated in genome-wide association studies (GWAS) cannot capture the effects of DNA variability on 'around the clock' glucoregulatory processes. Here we show that GWAS meta-analysis of glucose measurements under nonstandardized conditions (random glucose (RG)) in 476,326 individuals of diverse ancestries and without diabetes enables locus discovery and innovative pathophysiological observations. We discovered 120 RG loci represented by 150 distinct signals, including 13 with sex-dimorphic effects, two cross-ancestry and seven rare frequency signals. Of these, 44 loci are new for glycemic traits. Regulatory, glycosylation and metagenomic annotations highlight ileum and colon tissues, indicating an underappreciated role of the gastrointestinal tract in controlling blood glucose. Functional follow-up and molecular dynamics simulations of lower frequency coding variants in glucagon-like peptide-1 receptor (GLP1R), a type 2 diabetes treatment target, reveal that optimal selection of GLP-1R agonist therapy will benefit from tailored genetic stratification. We also provide evidence from Mendelian randomization that lung function is modulated by blood glucose and that pulmonary dysfunction is a diabetes complication. Our investigation yields new insights into the biology of glucose regulation, diabetes complications and pathways for treatment stratification.

Rare variants in the MECP2 gene in girls with central precocious puberty: a translational cohort study.

BACKGROUND: Identification of genetic causes of central precocious puberty have revealed epigenetic mechanisms as regulators of human pubertal timing. MECP2, an X-linked gene, encodes a chromatin-associated protein with a role in gene transcription. MECP2 loss-of-function mutations usually cause Rett syndrome, a severe neurodevelopmental disorder. Early pubertal development has been shown in several patients with Rett syndrome. The aim of this study was to explore whether MECP2 variants are associated with an idiopathic central precocious puberty phenotype. METHODS: In this translational cohort study, participants were recruited from seven tertiary centres from five countries (Brazil, Spain, France, the USA, and the UK). Patients with idiopathic central precocious puberty were investigated for rare potentially damaging variants in the MECP2 gene, to assess whether MECP2 might contribute to the cause of central precocious puberty. Inclusion criteria were the development of progressive pubertal signs (Tanner stage 2) before the age of 8 years in girls and 9 years in boys and basal or GnRH-stimulated LH pubertal concentrations. Exclusion criteria were the diagnosis of peripheral precocious puberty and the presence of any recognised cause of central precocious puberty (CNS lesions, known monogenic causes, genetic syndromes, or early exposure to sex steroids). All patients included were followed up at the outpatient clinics of participating academic centres. We used high-throughput sequencing in 133 patients and Sanger sequencing of MECP2 in an additional 271 patients. Hypothalamic expression of Mecp2 and colocalisation with GnRH neurons were determined in mice to show expression of Mecp2 in key nuclei related to pubertal timing regulation. FINDINGS: Between Jun 15, 2020, and Jun 15, 2022, 404 patients with idiopathic central precocious puberty (383 [95%] girls and 21 [5%] boys; 261 [65%] sporadic cases and 143 [35%] familial cases from 134 unrelated families) were enrolled and assessed. We identified three rare heterozygous likely damaging coding variants in MECP2 in five girls: a de novo missense variant (Arg97Cys) in two monozygotic twin sisters with central precocious puberty and microcephaly; a de novo missense variant (Ser176Arg) in one girl with sporadic central precocious puberty, obesity, and autism; and an insertion (Ala6_Ala8dup) in two unrelated girls with sporadic central precocious puberty. Additionally, we identified one rare heterozygous 3'UTR MECP2 insertion (36_37insT) in two unrelated girls with sporadic central precocious puberty. None of them manifested Rett syndrome. Mecp2 protein colocalised with GnRH expression in hypothalamic nuclei responsible for GnRH regulation in mice. INTERPRETATION: We identified rare MECP2 variants in girls with central precocious puberty, with or without mild neurodevelopmental abnormalities. MECP2 might have a role in the hypothalamic control of human pubertal timing, adding to the evidence of involvement of epigenetic and genetic mechanisms in this crucial biological process. FUNDING: Fundação de Amparo à Pesquisa do Estado de São Paulo, Conselho Nacional de Desenvolvimento Científico e Tecnológico, and the Wellcome Trust.

Missense3D-PPI: A Web Resource to Predict the Impact of Missense Variants at Protein Interfaces Using 3D Structural Data.

In 2019, we released Missense3D which identifies stereochemical features that are disrupted by a missense variant, such as introducing a buried charge. Missense3D analyses the effect of a missense variant on a single structure and thus may fail to identify as damaging surface variants disrupting a protein interface i.e., a protein-protein interaction (PPI) site. Here we present Missense3D-PPI designed to predict missense variants at PPI interfaces. Our development dataset comprised of 1,279 missense variants (pathogenic n = 733, benign n = 546) in 434 proteins and 545 experimental structures of PPI complexes. Benchmarking of Missense3D-PPI was performed after dividing the dataset in training (320 benign and 320 pathogenic variants) and testing (226 benign and 413 pathogenic). Structural features affecting PPI, such as disruption of interchain bonds and introduction of unbalanced charged interface residues, were analysed to assess the impact of the variant at PPI. The performance of Missense3D-PPI was superior to that of Missense3D: sensitivity 44 % versus 8% and accuracy 58% versus 40%, p = 4.23 × 10-16. However, the specificity of Missense3D-PPI was lower compared to Missense3D (84% versus 98%). On our dataset, Missense3D-PPI's accuracy was superior to BeAtMuSiC (p = 3.4 × 10-5), mCSM-PPI2 (p = 1.5 × 10-12) and MutaBind2 (p = 0.0025). Missense3D-PPI represents a valuable tool for predicting the structural effect of missense variants on biological protein networks and is available at the Missense3D web portal (http://missense3d.bc.ic.ac.uk).

Somatic mutations of CADM1 in aldosterone-producing adenomas and gap junction-dependent regulation of aldosterone production.

Aldosterone-producing adenomas (APAs) are the commonest curable cause of hypertension. Most have gain-of-function somatic mutations of ion channels or transporters. Herein we report the discovery, replication and phenotype of mutations in the neuronal cell adhesion gene CADM1. Independent whole exome sequencing of 40 and 81 APAs found intramembranous p.Val380Asp or p.Gly379Asp variants in two patients whose hypertension and periodic primary aldosteronism were cured by adrenalectomy. Replication identified two more APAs with each variant (total, n = 6). The most upregulated gene (10- to 25-fold) in human adrenocortical H295R cells transduced with the mutations (compared to wildtype) was CYP11B2 (aldosterone synthase), and biological rhythms were the most differentially expressed process. CADM1 knockdown or mutation inhibited gap junction (GJ)-permeable dye transfer. GJ blockade by Gap27 increased CYP11B2 similarly to CADM1 mutation. Human adrenal zona glomerulosa (ZG) expression of GJA1 (the main GJ protein) was patchy, and annular GJs (sequelae of GJ communication) were less prominent in CYP11B2-positive micronodules than adjacent ZG. Somatic mutations of CADM1 cause reversible hypertension and reveal a role for GJ communication in suppressing physiological aldosterone production.

Protein structure-based evaluation of missense variants: Resources, challenges and future directions.

We provide an overview of the methods that can be used for protein structure-based evaluation of missense variants. The algorithms can be broadly divided into those that calculate the difference in free energy (ΔΔG) between the wild type and variant structures and those that use structural features to predict the damaging effect of a variant without providing a ΔΔG. A wide range of machine learning approaches have been employed to develop those algorithms. We also discuss challenges and opportunities for variant interpretation in view of the recent breakthrough in three-dimensional structural modelling using deep learning.

The effect of adjusting LDL-cholesterol for Lp(a)-cholesterol on the diagnosis of familial hypercholesterolaemia.

BACKGROUND: Familial hypercholesterolaemia (FH) diagnostic tools help prioritise patients for genetic testing and include LDL-C estimates commonly calculated using the Friedewald equation. However, cholesterol contributions from lipoprotein(a) (Lp(a)) can overestimate 'true' LDL-C, leading to potentially inappropriate clinical FH diagnosis. OBJECTIVE: To assess how adjusting LDL-C for Lp(a)-cholesterol affects FH diagnoses using Simon Broome (SB) and Dutch Lipid Clinic Network (DLCN) criteria. METHODS: Adults referred to a tertiary lipid clinic in London, UK were included if they had undergone FH genetic testing based on SB or DLCN criteria. LDL-C was adjusted for Lp(a)-cholesterol using estimated cholesterol contents of 17.3%, 30% and 45%, and the effects of these adjustments on reclassification to 'unlikely' FH and diagnostic accuracy were determined. RESULTS: Depending on the estimated cholesterol content applied, LDL-C adjustment reclassified 8-23% and 6-17% of patients to 'unlikely' FH using SB and DLCN criteria, respectively. The highest reclassification rates were observed following 45% adjustment in mutation-negative patients with higher Lp(a) levels. This led to an improvement in diagnostic accuracy (46% to 57% with SB, and 32% to 44% with DLCN following 45% adjustment) through increased specificity. However all adjustment factors led to erroneous reclassification of mutation-positive patients to 'unlikely' FH. CONCLUSION: LDL-C adjustment for Lp(a)-cholesterol improves the accuracy of clinical FH diagnostic tools. Adopting this approach would reduce unnecessary genetic testing but also incorrectly reclassify mutation-positive patients. Health economic analysis is needed to balance the risks of over- and under-diagnosis before LDL-C adjustments for Lp(a) can be recommended.

GWYRE: A Resource for Mapping Variants onto Experimental and Modeled Structures of Human Protein Complexes.

Rapid progress in structural modeling of proteins and their interactions is powered by advances in knowledge-based methodologies along with better understanding of physical principles of protein structure and function. The pool of structural data for modeling of proteins and protein-protein complexes is constantly increasing due to the rapid growth of protein interaction databases and Protein Data Bank. The GWYRE (Genome Wide PhYRE) project capitalizes on these developments by advancing and applying new powerful modeling methodologies to structural modeling of protein-protein interactions and genetic variation. The methods integrate knowledge-based tertiary structure prediction using Phyre2 and quaternary structure prediction using template-based docking by a full-structure alignment protocol to generate models for binary complexes. The predictions are incorporated in a comprehensive public resource for structural characterization of the human interactome and the location of human genetic variants. The GWYRE resource facilitates better understanding of principles of protein interaction and structure/function relationships. The resource is available at http://www.gwyre.org.

A common TMPRSS2 variant has a protective effect against severe COVID-19.

BACKGROUND: The human protein transmembrane protease serine type 2 (TMPRSS2) plays a key role in SARS-CoV-2 infection, as it is required to activate the virus' spike protein, facilitating entry into target cells. We hypothesized that naturally-occurring TMPRSS2 human genetic variants affecting the structure and function of the TMPRSS2 protein may modulate the severity of SARS-CoV-2 infection. METHODS: We focused on the only common TMPRSS2 non-synonymous variant predicted to be damaging (rs12329760 C>T, p.V160M), which has a minor allele frequency ranging from 0.14 in Ashkenazi Jewish to 0.38 in East Asians. We analysed the association between the rs12329760 and COVID-19 severity in 2,244 critically ill patients with COVID-19 from 208 UK intensive care units recruited as part of the GenOMICC (Genetics Of Mortality In Critical Care) study. Logistic regression analyses were adjusted for sex, age and deprivation index. For in vitro studies, HEK293 cells were co-transfected with ACE2 and either TMPRSS2 wild type or mutant (TMPRSS2V160M). A SARS-CoV-2 pseudovirus entry assay was used to investigate the ability of TMPRSS2V160M to promote viral entry. RESULTS: We show that the T allele of rs12329760 is associated with a reduced likelihood of developing severe COVID-19 (OR 0.87, 95%CI:0.79-0.97, p = 0.01). This association was stronger in homozygous individuals when compared to the general population (OR 0.65, 95%CI:0.50-0.84, p = 1.3 × 10-3). We demonstrate in vitro that this variant, which causes the amino acid substitution valine to methionine, affects the catalytic activity of TMPRSS2 and is less able to support SARS-CoV-2 spike-mediated entry into cells. CONCLUSION: TMPRSS2 rs12329760 is a common variant associated with a significantly decreased risk of severe COVID-19. Further studies are needed to assess the expression of TMPRSS2 across different age groups. Moreover, our results identify TMPRSS2 as a promising drug target, with a potential role for camostat mesilate, a drug approved for the treatment of chronic pancreatitis and postoperative reflux esophagitis, in the treatment of COVID-19. Clinical trials are needed to confirm this.

Lipid lowering with inclisiran: a real-world single-centre experience.

OBJECTIVE: The reduction in circulating low-density lipoprotein cholesterol (LDL-c) is the primary aim of lipid-lowering therapies as a method of atherosclerotic cardiovascular disease risk reduction. Inclisiran is a new and potent lipid-lowering drug that is shown to be effective in reducing LDL-c in randomised controlled trials, however, real-world data of its use are not yet known. We sought to analyse the early effects of this drug in a tertiary centre lipid and cardiovascular risk clinic. METHODS: We performed a retrospective analysis of the first 80 patients who received a single dose of inclisiran at our lipid clinic between 1 December 2021 and 1 September 2022. Data were collected using electronic healthcare records. Baseline blood tests were taken prior to start of treatment and were repeated at 2 months follow-up. Data on adverse events were also recorded. RESULTS: At 2 months after treatment initiation, mean baseline LDL-c fell from 3.5±1.1 mmol/L by 48.6% to 1.8±1.0 mmol/L and total cholesterol from 5.7±1.3 mmol/L by 33.3% to 3.8±1.1 mmol/L (both p<0.0001). Mean high-density lipoprotein-c rose by 7.7% to 1.4±0.4 mmol/L (p=0.02) and median triglycerides fell by 31.3% to 1.1 mmol/L (IQR 0.9-2) (p=0.001). Adverse events (injection site reaction, fatigue and headache) were recorded in three patients and all had self-resolved by time of follow-up. CONCLUSION: Inclisiran use in line with National Institute for Health and Care Excellence guidelines led to significant lowering of LDL-c at 2 months, with efficacy similar to that reported in trials with good tolerability.

The AlphaFold Database of Protein Structures: A Biologist's Guide.

AlphaFold, the deep learning algorithm developed by DeepMind, recently released the three-dimensional models of the whole human proteome to the scientific community. Here we discuss the advantages, limitations and the still unsolved challenges of the AlphaFold models from the perspective of a biologist, who may not be an expert in structural biology.

Whole exome sequencing identifies deleterious rare variants in CCDC141 in familial self-limited delayed puberty.

Developmental abnormalities of the gonadotropin-releasing hormone (GnRH) neuronal network result in a range of conditions from idiopathic hypogonadotropic hypogonadism to self-limited delayed puberty. We aimed to discover important underlying regulators of self-limited delayed puberty through interrogation of GnRH pathways. Whole exome sequencing (WES) data consisting of 193 individuals, from 100 families with self-limited delayed puberty, was analysed using a virtual panel of genes related to GnRH development and function (n = 12). Five rare predicted deleterious variants in Coiled-Coil Domain Containing 141 (CCDC141) were identified in 21 individuals from 6 families (6% of the tested cohort). Homology modeling predicted all five variants to be deleterious. CCDC141 mutant proteins showed atypical subcellular localization associated with abnormal distribution of acetylated tubulin, and expression of mutants resulted in a significantly delayed cell migration, demonstrated in transfected HEK293 cells. These data identify mutations in CCDC141 as a frequent finding in patients with self-limited delayed puberty. The mis-localization of acetylated tubulin and reduced cell migration seen with mutant CCDC141 suggests a role of the CCDC141-microtubule axis in GnRH neuronal migration, with heterozygous defects potentially impacting the timing of puberty.

Missense3D-DB web catalogue: an atom-based analysis and repository of 4M human protein-coding genetic variants.

The interpretation of human genetic variation is one of the greatest challenges of modern genetics. New approaches are urgently needed to prioritize variants, especially those that are rare or lack a definitive clinical interpretation. We examined 10,136,597 human missense genetic variants from GnomAD, ClinVar and UniProt. We were able to perform large-scale atom-based mapping and phenotype interpretation of 3,960,015 of these variants onto 18,874 experimental and 84,818 in house predicted three-dimensional coordinates of the human proteome. We demonstrate that 14% of amino acid substitutions from the GnomAD database that could be structurally analysed are predicted to affect protein structure (n = 568,548, of which 566,439 rare or extremely rare) and may, therefore, have a yet unknown disease-causing effect. The same is true for 19.0% (n = 6266) of variants of unknown clinical significance or conflicting interpretation reported in the ClinVar database. The results of the structural analysis are available in the dedicated web catalogue Missense3D-DB ( http://missense3d.bc.ic.ac.uk/ ). For each of the 4 M variants, the results of the structural analysis are presented in a friendly concise format that can be included in clinical genetic reports. A detailed report of the structural analysis is also available for the non-experts in structural biology. Population frequency and predictions from SIFT and PolyPhen are included for a more comprehensive variant interpretation. This is the first large-scale atom-based structural interpretation of human genetic variation and offers geneticists and the biomedical community a new approach to genetic variant interpretation.

Ligand-Specific Factors Influencing GLP-1 Receptor Post-Endocytic Trafficking and Degradation in Pancreatic Beta Cells.

The glucagon-like peptide-1 receptor (GLP-1R) is an important regulator of blood glucose homeostasis. Ligand-specific differences in membrane trafficking of the GLP-1R influence its signalling properties and therapeutic potential in type 2 diabetes. Here, we have evaluated how different factors combine to control the post-endocytic trafficking of GLP-1R to recycling versus degradative pathways. Experiments were performed in primary islet cells, INS-1 832/3 clonal beta cells and HEK293 cells, using biorthogonal labelling of GLP-1R to determine its localisation and degradation after treatment with GLP-1, exendin-4 and several further GLP-1R agonist peptides. We also characterised the effect of a rare GLP1R coding variant, T149M, and the role of endosomal peptidase endothelin-converting enzyme-1 (ECE-1), in GLP1R trafficking. Our data reveal how treatment with GLP-1 versus exendin-4 is associated with preferential GLP-1R targeting towards a recycling pathway. GLP-1, but not exendin-4, is a substrate for ECE-1, and the resultant propensity to intra-endosomal degradation, in conjunction with differences in binding affinity, contributes to alterations in GLP-1R trafficking behaviours and degradation. The T149M GLP-1R variant shows reduced signalling and internalisation responses, which is likely to be due to disruption of the cytoplasmic region that couples to intracellular effectors. These observations provide insights into how ligand- and genotype-specific factors can influence GLP-1R trafficking.

LGR4 deficiency results in delayed puberty through impaired Wnt/ß-catenin signaling.

The initiation of puberty is driven by an upsurge in hypothalamic gonadotropin-releasing hormone (GnRH) secretion. In turn, GnRH secretion upsurge depends on the development of a complex GnRH neuroendocrine network during embryonic life. Although delayed puberty (DP) affects up to 2% of the population, is highly heritable, and is associated with adverse health outcomes, the genes underlying DP remain largely unknown. We aimed to discover regulators by whole-exome sequencing of 160 individuals of 67 multigenerational families in our large, accurately phenotyped DP cohort. LGR4 was the only gene remaining after analysis that was significantly enriched for potentially pathogenic, rare variants in 6 probands. Expression analysis identified specific Lgr4 expression at the site of GnRH neuron development. LGR4 mutant proteins showed impaired Wnt/ß-catenin signaling, owing to defective protein expression, trafficking, and degradation. Mice deficient in Lgr4 had significantly delayed onset of puberty and fewer GnRH neurons compared with WT, whereas lgr4 knockdown in zebrafish embryos prevented formation and migration of GnRH neurons. Further, genetic lineage tracing showed strong Lgr4-mediated Wnt/ß-catenin signaling pathway activation during GnRH neuron development. In conclusion, our results show that LGR4 deficiency impairs Wnt/ß-catenin signaling with observed defects in GnRH neuron development, resulting in a DP phenotype.

A polygenic biomarker to identify patients with severe hypercholesterolemia of polygenic origin.

BACKGROUND: Severe hypercholesterolemia (HC, LDL-C > 4.9 mmol/L) affects over 30 million people worldwide. In this study, we validated a new polygenic risk score (PRS) for LDL-C. METHODS: Summary statistics from the Global Lipid Genome Consortium and genotype data from two large populations were used. RESULTS: A 36-SNP PRS was generated using data for 2,197 white Americans. In a replication cohort of 4,787 Finns, the PRS was strongly associated with the LDL-C trait and explained 8% of its variability (p = 10-41 ). After risk categorization, the risk of having HC was higher in the high- versus low-risk group (RR = 4.17, p < 1 × 10-7 ). Compared to a 12-SNP LDL-C raising score (currently used in the United Kingdom), the PRS explained more LDL-C variability (8% vs. 6%). Among Finns with severe HC, 53% (66/124) versus 44% (55/124) were classified as high risk by the PRS and LDL-C raising score, respectively. Moreover, 54% of individuals with severe HC defined as low risk by the LDL-C raising score were reclassified to intermediate or high risk by the new PRS. CONCLUSION: The new PRS has a better predictive role in identifying HC of polygenic origin compared to the currently available method and can better stratify patients into diagnostic and therapeutic algorithms.

Identification of disease-associated loci using machine learning for genotype and network data integration.

MOTIVATION: Integration of different omics data could markedly help to identify biological signatures, understand the missing heritability of complex diseases and ultimately achieve personalized medicine. Standard regression models used in Genome-Wide Association Studies (GWAS) identify loci with a strong effect size, whereas GWAS meta-analyses are often needed to capture weak loci contributing to the missing heritability. Development of novel machine learning algorithms for merging genotype data with other omics data is highly needed as it could enhance the prioritization of weak loci. RESULTS: We developed cNMTF (corrected non-negative matrix tri-factorization), an integrative algorithm based on clustering techniques of biological data. This method assesses the inter-relatedness between genotypes, phenotypes, the damaging effect of the variants and gene networks in order to identify loci-trait associations. cNMTF was used to prioritize genes associated with lipid traits in two population cohorts. We replicated 129 genes reported in GWAS world-wide and provided evidence that supports 85% of our findings (226 out of 265 genes), including recent associations in literature (NLGN1), regulators of lipid metabolism (DAB1) and pleiotropic genes for lipid traits (CARM1). Moreover, cNMTF performed efficiently against strong population structures by accounting for the individuals' ancestry. As the method is flexible in the incorporation of diverse omics data sources, it can be easily adapted to the user's research needs. AVAILABILITY AND IMPLEMENTATION: An R package (cnmtf) is available at https://lgl15.github.io/cnmtf_web/index.html. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

PhyreRisk: A Dynamic Web Application to Bridge Genomics, Proteomics and 3D Structural Data to Guide Interpretation of Human Genetic Variants.

PhyreRisk is an open-access, publicly accessible web application for interactively bridging genomic, proteomic and structural data facilitating the mapping of human variants onto protein structures. A major advance over other tools for sequence-structure variant mapping is that PhyreRisk provides information on 20,214 human canonical proteins and an additional 22,271 alternative protein sequences (isoforms). Specifically, PhyreRisk provides structural coverage (partial or complete) for 70% (14,035 of 20,214 canonical proteins) of the human proteome, by storing 18,874 experimental structures and 84,818 pre-built models of canonical proteins and their isoforms generated using our in house Phyre2. PhyreRisk reports 55,732 experimentally, multi-validated protein interactions from IntAct and 24,260 experimental structures of protein complexes. Another major feature of PhyreRisk is that, rather than presenting a limited set of precomputed variant-structure mapping of known genetic variants, it allows the user to explore novel variants using, as input, genomic coordinates formats (Ensembl, VCF, reference SNP ID and HGVS notations) and Human Build GRCh37 and GRCh38. PhyreRisk also supports mapping variants using amino acid coordinates and searching for genes or proteins of interest. PhyreRisk is designed to empower researchers to translate genetic data into protein structural information, thereby providing a more comprehensive appreciation of the functional impact of variants. PhyreRisk is freely available at http://phyrerisk.bc.ic.ac.uk.

Can Predicted Protein 3D Structures Provide Reliable Insights into whether Missense Variants Are Disease Associated?

Knowledge of protein structure can be used to predict the phenotypic consequence of a missense variant. Since structural coverage of the human proteome can be roughly tripled to over 50% of the residues if homology-predicted structures are included in addition to experimentally determined coordinates, it is important to assess the reliability of using predicted models when analyzing missense variants. Accordingly, we assess whether a missense variant is structurally damaging by using experimental and predicted structures. We considered 606 experimental structures and show that 40% of the 1965 disease-associated missense variants analyzed have a structurally damaging change in the mutant structure. Only 11% of the 2134 neutral variants are structurally damaging. Importantly, similar results are obtained when 1052 structures predicted using Phyre2 algorithm were used, even when the model shares low (<40%) sequence identity to the template. Thus, structure-based analysis of the effects of missense variants can be effectively applied to homology models. Our in-house pipeline, Missense3D, for structurally assessing missense variants was made available at http://www.sbg.bio.ic.ac.uk/~missense3d.

EAP1 regulation of GnRH promoter activity is important for human pubertal timing.

The initiation of puberty is orchestrated by an augmentation of gonadotropin-releasing hormone (GnRH) secretion from a few thousand hypothalamic neurons. Recent findings have indicated that the neuroendocrine control of puberty may be regulated by a hierarchically organized network of transcriptional factors acting upstream of GnRH. These include enhanced at puberty 1 (EAP1), which contributes to the initiation of female puberty through transactivation of the GnRH promoter. However, no EAP1 mutations have been found in humans with disorders of pubertal timing. We performed whole-exome sequencing in 67 probands and 93 relatives from a large cohort of familial self-limited delayed puberty (DP). Variants were analyzed for rare, potentially pathogenic variants enriched in case versus controls and relevant to the biological control of puberty. We identified one in-frame deletion (Ala221del) and one rare missense variant (Asn770His) in EAP1 in two unrelated families; these variants were highly conserved and potentially pathogenic. Expression studies revealed Eap1 mRNA abundance in peri-pubertal mouse hypothalamus. EAP1 binding to the GnRH1 promoter increased in monkey hypothalamus at the onset of puberty as determined by chromatin immunoprecipitation. Using a luciferase reporter assay, EAP1 mutants showed a reduced ability to trans-activate the GnRH promoter compared to wild-type EAP1, due to reduced protein levels caused by the Ala221del mutation and subcellular mislocation caused by the Asn770His mutation, as revealed by western blot and immunofluorescence, respectively. In conclusion, we have identified the first EAP1 mutations leading to reduced GnRH transcriptional activity resulting in a phenotype of self-limited DP.